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1. Introduction 

In the last decades the phase diagram of quantum chromodynamics (QCD) has been a very ac- 
tive field of research. For small quark chemical potential, but non- vanishing temperature enormous 
progress has been achieved, which is mainly due to the fact that both, first principle continuum 
methods, see e.g. [1, 2, 3, 4, 5], as well as lattice QCD, see e.g. [6, 7, 8], have been applied very 
successfully. Nonetheless, even at vanishing density, some questions still remain open, in partic- 
ular in the vicinity of the chiral and confinement-deconfmement phase transitions. This concerns 
for example the trace anomaly in this region. Indeed, for temperatures around the deconfinement- 
confinement phase transition T <3T C continuum methods so far lack quantitative precision. This 
temperature regime is governed by fully non-perturbative physics, and requires the application of 
non-perturbative techniques. 

In continuum approaches to QCD, it is mainly the glue sector which requires the most effort, in 
particular conceptually but also technically. In the past two decades much progress has been made 
in order to understand Yang-Mills theory in the Landau gauge at vanishing temperature, leading to 
quantitative precision for correlation functions, see [9]. At finite temperature, however, apart from 
conceptual issues, the same quantitative precision is lacking. 

In the present contribution we report on work published in [10]. We aim at a quantitative 
computation of the Yang-Mills thermal propagators, for a review on thermal gluons see [11]. Here 
we utilise the functional renormalisation group (FRG), for QCD-related reviews see [1, 12, 13, 
14, 15, 16, 17, 18]. The FRG incorporates non-perturbative effects by a successive integration of 
fluctuations related to a given momentum scale, and hence is applicable to all temperatures. In the 
presented work we will focus on the pure Yang-Mills part of QCD, which has so far been studied in 
the framework of functional methods by the help of Dyson-Schwinger equations [11, 19, 20]. We 
compute fully dressed, non-perturbative ghost and gluon propagators, as well as the ghost-gluon 
vertex at finite temperature. The electric propagator shows the typical thermal screening above 
the confinement-deconfinement temperature T c in quantitative agreement with the lattice results, 
see [11, 20, 21, 22, 23, 24, 25]. It also shows an enhancement below T c which is qualitatively 
similar to the lattice findings. The magnetic propagator shows thermal scaling and tends towards 
the three-dimensional one for large temperatures. As for the electric propagator, it agrees well with 
the lattice. 

2. Yang-Mills Green Functions with the Functional Renormalisation Group 

The FRG is an exact method that is derived directly from the functional integral. It constitutes a 
renormalisation group equation [26] for the scale-dependent effective action r\, which interpolates 
between the classical action in the ultraviolet and the full effective action in the infrared by variation 
of an infrared cutoff, which suppresses fluctuations with momenta smaller than the infrared scale 
k. In the limit of k -> one is left with the full quantum effective action. 

The classical gauge fixed SU (N c ) Yang-Mills action in Euclidean space in the Landau gauge 
is given by 




(2.1) 
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with the ghost fields c a ,c a . The field strength tensor Ffi y and the covariant derivative D ab are defined 
in terms of the gluons A^, the coupling g and the structure constants f abc of the gauge group, 

F; v = d^A a v - d v A^+gf abc A b ^A c v , Df = 8 ab d^+gr b Al . (2.2) 

In the FRG-framework, infrared fluctuations below the scale k are suppressed by a modification of 
the propagator in the classical action of the form S — )► S + AS*, with 

AS *= / 7^4^<v(/)4+ / T^tf^V. (2.3) 
J (27T) z J (271) 

The regulator functions R ab v and R ab carry the same Lorentz and colour structure as the corre- 
sponding field propagators and are proportional to the dimensionless shape function r such that 
schematically R(p 2 ,k 2 ) ~ p 2 r(p 2 /k 2 ). Note that the ghost regulator is negative related to the un- 
physical dispersion of the ghost. In the following we use regulators of the form given in Fig. 1. The 
flow equation for the scale-dependent Yang-Mills effective action r\ at finite temperature is given 
by 

d t T k [<p] = ^TrG[<p](p,p) ■ d t R k {p) • (2-4) 

The trace stands for summation over species of fields with a relative minus sign for fermions, a 
summation over Lorentz and colour indices, and the integration/summation of the loop momenta. 
We also introduced the compact notation <p = (A,c,c). In eq. (2.4) G[(p] is the short notation 
for the fully dressed field dependent propagator. At imaginary time in the Matsubara formalism 
the integration measure at finite temperature turns the integration over the temporal momentum 
component into a sum over Matsubara frequencies, 

'(^4 - ^1/(03. ™* * = <") 

The functional flow equation for the effective action (2.4) can be illustrated in a diagrammatic form, 
which is given in Fig. 2. The minus sign in front of the ghost loops originates in the fermionic 
nature of the ghost. The flow equations for the propagators can be obtained by taking functional 
derivatives with respect to two gluons for the gluon propagator, and with respect to a ghost and 
an antighost for the ghost-propagator. These equations for the Yang-Mills two-point functions are 
depicted in Fig. 3. They only contain one-loop terms in full propagators and vertices. We emphasise 
that the flow equation depicted in Fig. 3 is exact, no higher loop terms are missing. This originates 
in the one-loop form of (2.4), see also Fig. 2. 



/< 
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3. Thermal Fluctuations 

The integration of the flow equation (2.4) involves both, quantum as well as thermal fluctu- 
ations. However, the treatment of the thermal fluctuations can be done separately via the help of 
purely thermal flows, which are constructed as the difference of the flows at finite and zero temper- 
ature, [13, 27], 



The quantity Ar r? ^[0] = r\[0]| r — r^[0]| r=o accounts for the difference between the effective 
action at finite and zero temperature. For the implementation of this procedure the vacuum physics 
is used as the starting point. A given set of correlation functions at vanishing temperature and 
scale k = is integrated by the help of the flow (2.4) up to a scale A, which is chosen such, 
that all thermal fluctuations are suppressed for the considered temperature T. Thus at this scale the 
correlation functions are not altered at leading order by switching on the temperature. With the flow 
at finite temperature we can flow down to k = again, this time integrating thermal fluctuations 
as well. Note that within this procedure only the difference d t ATx^ i- e - oni y thermal effects, are 
sensitive to the chosen approximation. 

The idea of the thermal flow is illustrated in a heuristic plot in Fig. 4. However, note that 
for increasing renormalisation group scale k the decay AT^[(/?i, ...,/?„) —> is exponential with 
exponent k/T but only polynomially in pl/k 2 , which is an important issue for numerical accuracy 
of thermodynamic quantities. This polynomial suppression is due to the fact, that only the loop- 
momenta are constrained by the regulator q < k, but the external momenta p are not limited. For 
large momenta p » k the flow factorises at leading order, see Fig. 5, which weakens the locality 
of the flow. However, the locality of a flow can be restored by a ^-dependent reparametrisation 
of the fields, which can be done in such a way, that the flow of the two-point functions is strictly 
zero for external momenta larger than a certain scale p > Xk, [15, 10]. Note that the contributions 
stemming from this region are not dropped, but rather stored in a reparametrisation function of the 
fields. This reparametrisation must be undone at the end of the computation. In the following we 
will drop any reference to the reparametrisation of fields and refer the interested reader to [10]. 

4. Approximation 

The flow of an arbitrary n-point function depends itself on m-point functions with m < n + 2, 
see e.g. that of the propagators depicted in Fig. 3. As a consequence, the flow in eq. (2.4) gener- 
ates an infinite hierarchy of coupled integro-differential equations. Within numerical applications 



Figure 2: Functional flow for the effective action. Lines with filled circles denote fully dressed field depen- 
dent propagators G[(p]. Crossed circles denote the regulator insertion d t Rk- 



dATrM = l^G[(j)} ■ d t R\ T - iTrG[0] • d t R 



T=0 ' 



(3.1) 
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Figure 3: Flow equations for the ghost and gluon propagators. Vertices with filled circles denote fully 
dressed vertices. All propagators are fully dressed, the filled circles for the internal ones have been dropped 
for the sake of clarity. Crossed circles denote the regulator insertion d t Rk. 




Figure 4: Flow at vanishing temperature from k = to k — A, and flow at T ^ from k = A to k = 
with r/Acl. The flow is described in theory space, and the axes label (orthogonal) couplings/observables 
Gi which serve as expansion coefficients of the effective action, e.g. G\ = T^ 2 \p = 0). The flows start to 
deviate at k « 2kT. 



this system must be closed by projecting it on a finite set of flow equations while keeping the rel- 
evant physics. In other words, the contributing correlation functions must already comprise the 
relevant physics. For the propagators we keep the full propagators and work with self-consistent 
approximations to the vertices which respect the renormalisation group properties of the vertices. 
Furthermore, we use the flow equation for the ghost-gluon vertex at the symmetric point with 

P 2 =k 2 . 

In the Landau gauge the gluon is a purely transversal vector particle. At finite temperature 
the heat bath singles out a preferred rest frame defined by a time-like vector n^. This leaves us 
with one (3d) longitudinal and one transversal tensor structure, both being transverse in Ad. The 
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p + q P 



Figure 5: Factorisation in leading order for large momenta for the first diagram (cutted gluon line) in the 
flow of the (inverse) ghost propagator given in Fig. 3. The triangle stands for the product of the two vertices 
at q = and reads —p\xP^f acd f bc d • 



corresponding projection operators are given by 

U Iv(p) = Siiv-PiiPv/p 2 , 

p ?iv(P0iP) = (l-8 ll0 )(l-8 v0 )(8 llv -p ll p v /p 2 ) , 

Ptv(po,P) = nj v (p) -p£ v (po,p). (4.1) 

At vanishing temperature the gluon can be parametrised with one scalar function. At non- vanishing 
temperature the wave-function renormalisations for the longitudinal propagator Z L and transverse 
propagator Z T differ and have to be taken into account separately. The ghost has no Lorenz struc- 
ture and can be described fully in terms of one scalar function also at non- vanishing temperature. 
Leaving out the identity in colour-space the parametrisation is done according to 

( r S) Aiv (/ ? o,p) = z l ( POi p)p 2 pI; v (po,p), 

( r S) Aiv (/ ? o,p) = Z T (p ,p)p 2 P^ v (p () ,p), 

T?\p Q ,p)=Z c {p Q ,p)p\ (4.2) 

where we emphasise that the wave-function renormalisations depend on the temporal momentum 
component po and the spatial ones p separately. 

The flow equations for the two-point functions depend on the two-, three- and four-point func- 
tions. In particular we have tadpole diagrams which depend on the ghost-ghost and ghost-gluon 
scattering vertices r~ c ^ c and Ti^. Albeit these vertices are absent on the classical level they 
become relevant in the non-perturbative regime. They can be taken into account via a DSE resum- 
mation of the vertices in the flow. This turns the flow into the total scale derivative of the DSE, 
see Fig. 7(a). This reflects the fact that, in general, the flow equation is the differential form of the 
Dyson-Schwinger equation. We remark that this usually would introduce the necessity of renor- 
malising the corresponding DSE. In the present parameterisation this is avoided due to locality, 
for details see [10]. As a result the total scale derivative of the DSE for the ghost two-point func- 
tion in the presence of a regulator term is equivalent to the flow equation, but is free of four-point 
functions, for details see [10]. The scale derivative acting on the dressed propagator gives 

d t G [0] = -G [0] • d t (rP) [<j>] + V) • G [0] . (4.3) 
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x? 

Figure 6: Ghost-gluon vertex. 



For the vertices we introduce a parametrisation that naturally captures the renormalisation group 
behaviour, for details see [28]. With the wave-function renormalisation of the ghost and gluon we 
write 

r^(p u ...,p n ) = flzf(pi) n P u..,Pn) ■ (4.4) 

i=l 

The Z ( j > (p) -factors carry the RG- scaling of the vertex as well as the momentum dependence, thus 
they take into account potential kinematic singularities. The 2F are renormalisation group invariant 
tensors, which carry the canonical momentum dimension as well as the tensor and colour structure. 
For the ghost propagator we have sketched the idea how the ghost-ghost and ghost-gluon scattering 
vertices can be absorbed via a DSE resummation. With a similar resummation in the flow of the 
ghost-gluon vertex in Fig. 6 we can fully resolve its RG running at the symmetric point p 2 = k 2 . 
We use the parametrisation 

Z&Al'P) = ^Ac-[sil(q,p)]f c = z k ,, Ac iq^f abc , (4.5) 

§ 

where g is the classical coupling. The RG invariant factor Zk,cAc defines the running coupling 
a s = zlc Ac / (47r). Also for the three- and four-gluon vertices we restrict ourselves to the classical 
tensor and colour structure. The parametrisation is done according to eq. (4.4), with 

^ A 3 = z M 3 -5$ , = z M 4 • (4.6) 

For large momenta the couplings z^ A 3, z kA A relate to the ghost-gluon dressing z^cAc I n tne in- 
frared the couplings are suppressed strongly. This is taken into account via functions that approach 
Z3 = Z4 = 1 for k » Aqcd, but suppress the gluonic vertices in the infrared. The couplings are given 
by z^a 3 = Z3Zk,cAa Zk^A 4 = z 4 z lcAc ^ n Edition to that we regularise the vertex prefactors in such 
a way that we freeze the Za for scales smaller than the minimal turning point in the wave-function 
renormalisation. The naive ansatz of the pure wave-function renormalisation would yield unphys- 
ical IR divergent flows, which would spoil the IR suppression of the purely gluonic vertices. The 
above approximation scheme is depicted in Fig. 7, including the DSE resummations. 

5. Results for Propagators and the Ghost-Gluon Vertex 

In this section we present the results for the ghost and gluon propagators and the ghost-gluon 
vertex. The temperature is given in lattice units. Only below the typical temperature scale of 2kT 
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(a) Flows for the two-point functions in the given truncation. 




(b) Flow for the ghost-gluon vertex in the given truncation. 

Figure 7: Yang-Mills flows for propagators (a) and ghost-gluon vertex (b) the approximation discussed 
above. The flows for the ghost propagator and the ghost-gluon vertex are DSE-resummed and the ghost- 
tadpole in the gluon equation is neglected. 

we have significant thermal effects on the momentum dependence. In turn, above this scale the 
temperature fluctuations are suppressed and all propagators tend towards their vacuum counterparts 
at vanishing temperature. This also holds true for the ghost-gluon vertex. 

The most significant effect can be seen for the chromoelectric and chromomagnetic gluon 
propagators, the components of the propagator longitudinal and transversal to the heat bath re- 
spectively. The zero mode of the longitudinal gluon propagator at various temperatures is given in 
Fig. 8(a) as a function of spatial momentum. 

For low temperatures T < 150 MeV we see an enhancement of the longitudinal propagator. 
Such an enhancement is also seen on the lattice, [1 1, 20, 21, 22, 23, 24, 25]. For higher temperatures 
the longitudinal propagator is suppressed relative to the gluon propagator at vanishing temperature. 
This is the expected behaviour caused by the Debye screening mass due to the thermal screening 
of the chromoelectric gluon. For asymptotically high temperatures » T c the chromoelectric gluon 
decouples. The onset of this behaviour at about T « 150 MeV is earlier as in the respective lattice 
computations [11, 20, 21, 22, 23, 24, 25], where the thermal decoupling takes place for tempera- 
tures larger than the critical temperature. 

In order to quantitatively capture this behaviour we have to extend our present truncation with 
a self-consistent inclusion of the Polyakov loop background as well as a better resolution of the 
purely gluonic vertices for momenta and frequencies below Aqcd T c . Above the confinement- 
deconfinement scale we see quantitative agreement with the lattice results. 

The transversal mode is not enhanced for small temperatures in clear distinction to the longi- 
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(c) Zero mode of the ghost wave function renormalisation. 



(d) Zero mode of the ghost propagator 



4.3 
4.2 




fif. 
if/ 

i. • / 

f 



z cAc 



T = MeV 
T = 100 MeV 
T = 200 MeV 
T = 300 MeV 
T = 500 MeV 
T = 1 GeV 



0.5 



1.5 2 2.5 

k [GeV] 



3.5 



(e) Ghost-gluon vertex coupling. 
Figure 8: Yang-Mills propagators and the ghost-gluon vertex at various temperatures. 



tudinal mode. It is monotonously decreased with temperature, see Fig. 8(a). Moreover, it develops 
a clear peak at about 500 MeV. This can be linked to positivity violation which has to be present 
for the transversal mode as in the high temperature limit it describes the remaining dynamical glu- 
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ons of three-dimensional Yang-Mills theory in the Landau gauge. The infrared bending is more 
pronounced as that of respective lattice results, which already holds at vanishing temperature. It 
is another choice for the decoupling solution as discussed in [9]. Moreover, its strength may also 
depend on the quantitative precision achieved for the gluonic vertices at these scales. In turn, for 
larger momenta the transversal propagator agrees well with the respective lattice propagator. 

The ghost shows only a small temperature dependence in contradistinction to the gluonic prop- 
agators. The temperature dependence is hardly evident on the level of the propagator Fig. 8(d), but 
can be resolved on the level of the wave-function renormalisation Fig. 8(c). The wave-function 
renormalisation is slightly suppressed, which corresponds to a successive enhancement of the 
ghost propagator at finite temperature. The enhancement of the ghost propagator is potentially 
self-amplifying as it feeds back into the flow of the ghost two-point function, see Fig. 7. This 
would cause a pole in the ghost propagator at some temperature if the ghost-gluon vertex would 
not change. However, in this case the flow of the latter is dominated by the ghost diagram, see 
Fig. 6. This non-trivial interplay of the flow of the ghost propagator with the flow of the ghost- 
gluon vertex leads to a self-stabilising system and prevents a further enhancement of the ghost. 
This effect is crucial for the stability of the solution of the Yang-Mills system at finite temperature. 
Indeed, for a constant vertex no solution could be obtained for intermediate temperatures T « T c 
and above. Thus we conclude that any reliable truncation must comprise direct thermal effects also 
in /i-point functions with n < 3. 

The self-stabilising property of the Yang-Mills system explained above is clearly seen in the 
temperature-dependence of the ghost-gluon vertex. The vertex is suppressed successively with the 
temperature, see Fig. 8(e), which in turn ensures the relatively mild change of the ghost propagator. 
Especially the sharp drop-off of the vertex at small scales k accounts for the smallness of the thermal 
fluctuations to the ghost propagator. 

In the following, we compare the propagators above with lattice results [21]. For this purpose, 
we scale the lattice data such that the lattice propagators at vanishing temperature match our nor- 
malisation at momenta p>l GeV. Take notice, that we did not use the lattice propagator as the 
initial condition, thus the deep infrared of the data deviates from our propagator already at zero 
temperature, which persists also in the propagators at finite temperature. Apart from that, there is 
quantitative agreement with the lattice data with respect to the (temperature dependent) momentum 
region, where the thermal effects appear. In Fig. 9(b) the transversal propagators are compared. 

Clearly, we match the lattice propagator, except for the strong bending of the FRG propagator 
in the infrared region. As already discussed above, this difference is a direct consequence of the 
deviation of ghost and gluon propagators at vanishing temperature in the deep infrared, and the 
ansatz for the gluonic vertices. However, the behaviour of the magnetic gluon agrees with the 
lattice observations for momenta above 500 MeV and all temperatures. In contrast to this, the 
electric gluon, see Fig. 9(a) on the lattice shows a qualitatively different behaviour for temperatures 
below and around the phase transition. Note also that although the longitudinal propagators agree 
for T = 0.361 T c w lOOMeV, it is exactly this region where the uncertainty due to the truncation for 
the gluonic vertices is potentially large, for a detailed discussion see [10]. Being aware of a possible 
truncation dependence in the deep infrared of the longitudinal propagator at low temperatures, we 
note that in the present truncation the electric gluon qualitatively shows the enhancement found on 
the lattice. Increasing the temperature this feature disappears, and we see a qualitatively different 
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(a) Comparison of the longitudinal propagator from above (b) Comparison of the transversal propagator from above 
with the corresponding lattice result. with the corresponding lattice result. 

Figure 9: Comparison of the gluon propagator with the lattice results 

effect for temperatures below T c . While the continuum result shows a strictly monotonic decreasing 
propagator, the counterpart on the lattice is enhanced in the confining regime, but reflects the phase 
transition in form of a rapid decrease at T c . Nevertheless, this deflection is expected to be missed 
in the present truncation, as the full Polyakov loop potential V(Ao), see [29, 30], is pivotal for the 
critical behaviour around the phase transition, see [24]. In a full calculation the inverse propagator 

(2) 

is proportional to the second derivative of the Polyakov loop F\ J L ~ V"(Ao), however in the work 
presented here this term was dropped. In any case this term is absent in the magnetic modes. 

6. Summary and Outlook 

We have computed temperature dependent Yang-Mills propagators in Landau gauge in the 
framework of the functional renormalisation group for T < 3T C . For this purpose we employ ther- 
mal flows, which encode thermal fluctuations in the difference of the full renormalisation group 
flow at non- vanishing and zero temperature. 

The chromoelectric propagator shows the expected Debye-screening for T > T c in quantitative 
agreement with the lattice results. For small temperatures it shows qualitatively the enhancement 
also seen on the lattice [11, 20, 21, 22, 23, 24, 25]. However, the significance of the lattice results 
so far as well as quantitative details are not settled yet. This concerns in particular the behaviour of 
the chromoelectric propagator at criticality, see [24]. We hope to add something to the clarification 
of this issue within an extension of the present work. There we resolve the dependence on the 
Polyakov loop as well as have a better grip on the gluonic vertices. The chromomagnetic propagator 
shows the expected thermal scaling and tends towards the three-dimensional gluon propagator in 
quantitative agreement with the lattice. 

The ghost propagator only shows a mild enhancement with temperature in agreement with the 
lattice. In contradistinction we see a strong thermal infrared suppression of the ghost-gluon vertex 
which increases with temperature. 
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At present we improve on the approximations and compute thermodynamic observables, e.g. 
the pressure. Furthermore, we extend our work to full QCD with 2 and 2+1 flavours in an extension 
of [31] also at finite density. 
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